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Abstract 

We propose a coarse-grained model for polymer chains and polymer networks based on the meso- 
scale dynamics. The model takes the internal degrees of freedom of the constituent polymer chains 
into account using memory functions and colored noises. We apply our model to dilute polymer 
solutions and polymer networks. A numerical simulation on a dilute polymer solution demonstrates 
the validity of the assumptions on the dynamics of our model. By applying this model to polymer 
networks, we find a transition in the dynamical behavior from an isolated chain state to a network 
state. 
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lierarchical structures of polymer gels are also utilized to develop many innovative 
3Q- 



A wide class of systems possesses hierarchical structures over large length scales. Typ- 
ical examples are critical liquids and soft matter such as polymers, surfactants, colloidal 
suspensions, and polymer gels. Due to the coexistence of the different length scales of the 
internal freedom degree, soft matter shows various anomorous and interesting phenomena 
including shear induced phase separation of polymer solutions jl| 2], and viscoelastic phase 
separations of entangled polymers jsj. In such point of view, the system of polymer gels 
is a somewhat interesting target, because it has widely distributed length scales, i.e. the 
size of monomers, the size of networks, and the size of the whose elastic body. Due to such 
hierarchical structures, many interesting phenomena such as swelling behavior coupled with 
inhomogeneity QjQ] and anomalous relaxation process in the dynamics of networks^] ?! oc- 
cur. The 
materials 

To realize these phenomena and to design these materials, it is important to model meso- 
scale structures and to bridge between these mesoscale structures and those on the larger 
and smaller scales. There have been many models for the polymer network systems on 
different length scales[ll]. However, due to the complexity of the polymer networks, these 
models are far from realistic especially on the intermediate length and time scales, which 
are important in understanding the experimental results and in materials designing. 

In this paper, we propose a coarse-grained model of polymer networks based on the meso- 
scale structures and dynamics. We perform numerical simulations of this model to verify 
the validity of our assumptions and modeling. 

Let us discuss our coarse-grained model for polymer networks. To reduce the degrees 
of freedom, we derive a set of dynamical equations for the polymer networks which are 
described in terms of the degrees of freedom of the crosslinkers. In Fig{TJ we show a schematic 
illustration of our reduction process of the degrees of freedom. Here we model the polymer 
network as a set of linear polymer chains connected by crosslinkers as is shown in Fig{T](a). 
We can describe this network topology using so-called adjacency matrix. To obtain explicit 
expressions for the memory functions and the colored noises based on the meso-scale polymer 
network model, we decompose polymer network into a set of linear polymer chains which 
are described by the Zimm model for the dilute polymer solution under hydrodynamic 
interaction(FigJTIb)) 12J. By a reduction of the degree of freedom of the dynamic equation 
for the linear polymer chain, we gain the memory functions and colored noises for the 
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polymer chain(Fig{T](c)). Then, we reconstruct the polymer network by connecting these 
chains using the adjacency matrix. With this procedure, we can express the motion of the 
crosslinkers without using information on the motion of the monomers (Figfljd)). After this 
reduction, the degrees of freedom of the monomers are reduced to the memory functions and 
colored noises [3]. 

Now, we discuss the dynamics of an isolated linear polymer chain, which is composed 
of N monomers. Let R(n, t), (n = 1, ...N — 1) be the position of i-ih monomer at time t, 
and define R+(i) = R(0, t), R-(t) = H(N,t), v± = d~R±/dt, and f-t as the positions and 
the velocities of crosslinkers at the chain ends and the forces acting on them. In this dilute 



polymer case, the dynamics of monomers are well described by the Zimm model[12|]. In the 
Zimm model, the equation of motion of individual monomer is expressed as 

^f = k L dmh ( n - m )^f + 1 ln + fn (1) 

(l»(*)i£(0> = 2k B Th(n - m)I5(t - t>) (2) 

where h(n — m) and t] n are mobility of the monomer and the noise acting on the n-th 
monomer and f n is defined as f n = i + NS(n) + f_N5(n — N). The mobility of the monomers 
is related to the noise r\ n by the fluctuation dissipation theorem eq.pj, where I is the 
unit tensor. In order to solve eq.([T]) formally, we introduce the Fourier series expansion of 
any vector variable Q n such as R n as Q n = Qo + 2 52 p =i Qp cos (^0 together with the 

pnn 

~w 



expressions for the Fourier coefficients Qo = ^ J^dnQ n and Q p = J dnQ n cos(— 




(a) (d) 
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FIG. 1: Schematic diagram of the reduction process of the degrees of freedom of polymer network, 
(a) A polymer network before the reduction, (b) A tagged single polymer chain in the polymer 
network, (c) The polymer chain described by memory functions ^Z 1 * 1 . (d) Polymer network after 
the reduction. 
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Here we neglected the sine modes in the Fourier series, because we focus on the dynamics 
of the end points of the polymer chain which do not excite the cosine modes. Then the 
equation of motion eqflTJ) is rewritten as 

pL p = (-\-k p X p + f p ) + f, p , (3) 

where k p = 6n 2 kBT(Nb 2 )~ 1 p 2 and ( p = (12ir 3 Nb 2 p) 1 ' 2 r] s . Here, the longest relaxation time 
t r — Ci/^i is called Rouse time. If we only focus on the dynamics of crosslinkers, the 
independent variables we need are R+(£) = R(0,t) and R_(£) = R(iV, t). Since eq.flBJ) 
is linear equation, we can directly integrate eq.([3]) over t. Then, we can sum up Fourier 
coefficients X p to derive R-t. As a result, the velocities of crosslinkers v± = dH±/dt are 
given by 

v± = J ds{g + (t - s)f ± (s) +g-(t- s)f T (s)} + £ ± , (4) 
where g±(t) and £±(£) are the memory functions and colored noises that are defined by 



9± 



(±i)V /2 / *p 3/ \' 

exp ( ) 



(5) 



Co ~ C P L Tr t r 

(i±(t)il(t')) = k B Tg + (t-t'), (6) 
(L(t)ZUt')) = k B Tg4t-t'). (7) 

Using eq.(jl]), the force acting on the crosslinkers which is caused by the motion of internal 
degrees of freedom of the polymer chain can be described as 

f± = J ds{$ + (t - s)v±(s) + $"(*- s)v T (s)} + £±. (8) 

The memory functions $ ± are given by 

$ ± {u) = ±g ± {u)/{gl{u)-g 2 _{u)), (9) 

where $ ± (t) are calculated from eqs.(!51) and numerically. In Fig {5J we show the result 
of the numerical evaluation of $ ± (t). We confirm that the memory functions $ ± (t) can be 
very well approximated by 

r R ) = TA ± (t/t )- a ± exp(-t/r R ), (10) 
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FIG. 2: Time dependence of the normalized memory kernels which are defined by $± = 
&±/A±tQ ± on (a) semi-log plot and on (b) log-log plot. The top (solid) curve is <& + , the middle 
(broken) curve is &~ , and the bottom (dotted) curve shows the single exponential function for 
comparison. 

where to are the characteristic time scales of the microscopic dynamics, and the numerical 
values of exponent a± are found to be a + = 0.66 and a_ = 0.15, respectively. 

Now, we derive a coarse-grained model for polymer networks based on the Langevin 
equation for a single chain eqs.fllj) (JSj) and ffTUj) . Let Rj(t) and Vj(t) be the position and 
the velocity of i-th crooslinker at time t, and Fi{t) be the force acting on it. We define 
iy the force acting on the i-th crosslinker imposed by the j-th crosslink through the chain 
connecting them. Since the system we discuss is composed of the linear polymer chains which 
obey the linear integral equation eq.([ED, we can derive the force Fj as Fj = Yl!j=o Aij^ij using 
linear superposition with the adjacency matrix Ay- for the polymer network. If i and j is 
connected, element of the adjacency matrix is defined as Ay = 1. On the other hand, 
Ay — 0, if i and j are unconnected. The equations of motion of the crosslinkers are given 
by 

|H, = v„ (11) 

dt 3 J -°° J 

M 1j =( y EA 4k *+-c)6y + Ai d *q. (13) 
^ k ' 

The memory functions $^ are defined as $^ = TA±(t/t )~ a± exp(— t/ry), and Ty is the 
maximum relaxation time of the chain between i-th and j-th crosslinkers. ( and k are 
friction coefficient of a crosslinker and spring constant of a chain. Because the viscosity 



term in eq. (fT2|) is in general much larger than the inertia term in polymer solution, we can 
neglect the latter term. Finally, we discretize eqs. ffTTj) and f|T2|) so that we can integrate 
them numerically by Euler method as 

Ri(t + At) = Vi(t)At + Ri(t) + 0(At 2 ), 



(14) 



v<(*) = J2 / ds M(s) 

■ oo 

^2 AtM k j{nAt)v 5 {t - nAt) + kVL kj R,j(t) + $ k (t) 



n=l 



{ZmJ(t + nAt)) 



2k B TI 



At 



l -1 



dsM,(s) 



5 0A t + M-MnAt) 



(15) 



(16) 



where the memory kernel Aiij(t) is related to the colored noise $,i{t) by the fluctuation 
dissipation theorem eq. (1161) . 

We calculate the mean-square displacement C(t), intermediate scattering function of the 
position of the crosslinkers g c (q, t), and the total scattering from both the crosslinkers and 
the monomers g(q, t) by molecular dynamics (MD) simulations using eqs. (H4l) and ffT5l) . 
Here, C(t), g c (q,t) and g(q, t) are defined by 

c(t) = ^(|R(t)-R(o)| 2 >, ( 17 ) 



g c (q,t) = ^(exp [iq • (R,-(t) - R^O))] ), 



g(q,t) = ^(^(q)^(q))(exp [iq • (Xy(t) - X„(0))]>, 

i,j,k,l 



(18) 
(19) 



where R» and Xy = (Rj + Rj)/2 are the position of the i-th crosslinker andthe position of 
center of mass of the chain between crosslinkers i and j. Here, 5V, (q) expresses the scattering 
from the chain connecting i-th and j-th crosslinkers, and defined as «SV,(q) = J2 m ex P(^ ( l" r m)) 
with r m being the position of m-th monomer that belongs to the chain between i-th and 
j-th crosslinkers. In eq. f|T9|) . we neglect correlation between the internal degrees of freedom 
of each chain and that of its center of mass. Then, (S'jj(q)S'^(q)) = (\Sij(q)\ 2 )5ikSji is 
approximated as Debye function, because the relaxation of (SV^q^^q)) is much faster 
than that of the correlation between the positions of centers of mass. 

To confirm the validity of our model, we first calculate C(t) for a single linear polymer 
chain, whose ends are labeled as 1 and 2. In this case, we set At = 0.001, A±tQ ± = 1.1, 




FIG. 3: (a) Auto-correlation function of the position of the one end point of a linear polymer chain. 
(b)Dynamic structure factor of a strand composed of a linear polymer chain. The solid line is 
result of coarse-grained model, and dotted line shows spring-beads model's one, where the dotted 
line is scrolled up for comparison. 

( = 0.25, k = 0.075, ksT = 1.0 and r 12 = 1.0, where we measure length and time in the 
units of x = R g and t — T\i — tr, where R g is the radius of gyration of the polymer chain. 
For this linear polymer, the adjacency matrix and the interaction matrix are given by 

The simulation results are shown in Figj3](a) along with the data of the non coarse-grained 
bead-spring model that will be described later. It is evident that in the short time regime 
(t <C t) the auto-correlation function C(t) behaves as C(t) ~ i 2 / 3 , while in the long time 
regime (t 3> i) this function can be approximated as C(t) ~ t. These results are consistent 

n n 

with the experimental study [14] and the scaling argument [151]. Therefore we confirm the 
validity of our model for a single polymer chain. The intermediate scattering function 
eq. (fl9|) for a dilute polymer solution is also calculated and is shown in Figj3](b). Here, we 
set q as \q\R g = 2tt . In the short time regime, the relaxation of #(q, t) is well described 
as the stretched exponential function g(q,t) ~ g(q, 0) exp(— Tit 2 / 3 ). On the other hand, in 
the long time regime, we can fit g(q, t) as g(q, t) ~ g(q, 0) exp(— r 2 t). From this numerical 
result, it is clear that the intermediate scattering function g(q, t) can be described only using 
the degrees of freedom of the crosslinkers. In this simulation, we also compare our model 
to the spring-bead model (non coarse-grained model). In Figj3l it is shown that our model 
correctly reproduces the non coarse-grained one, in spite of its degrees of freedom is smaller 
than the non coarse-graind onejlj 
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FIG. 4: Time dependence of intermediate scattering function of the position of the crosslinkers 
g c (t) in (a) short time regime and in (b) intermediate time regime. The solid line is the simulation 
result of polymer network, and the dotted line shows that of single polymer chain for comparison. 

Next, we apply this model to a polymer network to understand the hierarchical structure. 
In the simulation, we use the same parameters as those used in Fig|3j and we set the 
number of crosslinkers as iV = 1000. To describe connectivity of the polymer network 
systems, we prepare a set of N linear polymers, and connect any pair of the end points 

ity p = 0.0005, which corresponds 



15j. The interaction matrix flij is 



randomly with a probability p. Here, we set the probabi 
to the percolation point of this random network system 
given by flij = —(^2 k A ik )5ij + Aij. In FigHJ we show the result of simulations on g c (q, t) 
with |q|-R fl = 2tt. This numerical result shows that in short time regime (t/t < 0.1), the 
correlation function of the network is the same as the one for a single chain. On the other 
hand, in the intermediate time regime (0.1 < t/t < 1), the relaxation of this function is 
much slower than that of a single chain, and this is roughly approximated by a power low 
function as t~ a6 . This power law relaxation process can be explained by a simple scaling 
approach using the cluster distribution function P(n) ~ n a and the maximum relaxation 
time r max (n) ~ rfi at a percolation threshold. Since this percolated system can be classified 
as Bathe grid system, the cluster distribution function P{n) and the cluster radius R s is 
described in general as P(n) ~ n~ 25 and R s ~ n 0,25 . The relation between r max , R s 
and number of monomers form end to end of the cluster N are r max ~ iV 2 ~ Rf, and 
hence (3 = 1. In the present situation, because of finite size effect and hydrodynamic 
interaction between crosslinkers, the exponent (3 is slightly smaller than this expected value 
and is estimated as (3 ~ 5/6. Then, in the intermediate time regime g c (t) is described as 
9c(t) ~ / dn nP(ri)e~*/ Tmax( ' n ) ~ t( a + 2 )// 3 = t~ os , which reproduces the power low's behavior 
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shown in FigH^b). Thus, we conclude that the origin of the power law relaxation in the 
intermediate time scale (FigJH(b)) is the motion of percolated clusters. Similar phenomena 
were observed experimentally under the gelatin process[6j. 

In conclusion, we have proposed a coarse-grained model for polymer networks and showed 
its numerical results. The main results of the present work are summarized below, (i) 
We have derived a set of dynamic equations of polymer networks only using the degrees 
of freedom of the crosslinkers, and obtained the expressions of the quite general memory 
kernels and the random noises, (ii) The auto-correlation function C(t) of the crosslinkers at 
the both sides of a single polymer chain and their intermediate scattering function g(q, t) in 
dilute polymer solution were calculated numerically. The simulation results were consistent 
with the experimental and theoretical results [14J . (iii) We applied this model to a polymer 
network, and succeeded in smoothly connecting the dynamics of crosslinkers between "the 
internal motion of a single chain" and "the motion of percolated clusters". 
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